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ABSTRACT 

We present the first Hinode/EIS observations of 5 min quasi-periodic oscillations detected in a 
transition-region line (He II) and five coronal lines (Fe X, Fe XII, Fe XIII, Fe XIV, and Fe XV) at the 
footpoint of a coronal loop. The oscillations exist throughout the whole observation, characterized 
by a series of wave packets with nearly constant period, typically persisting for 4-6 cycles with a 
lifetime of 20-30 min. There is an approximate in-phase relation between Doppler shift and intensity 
oscillations. This provides evidence for slow magnetoacoustic waves propagating upwards from the 
transition region into the corona. We find that the oscillations detected in the five coronal lines are 
highly correlated, and the amplitude decreases with increasing temperature. The amplitude of Doppler 
shift oscillations decrease by a factor of about 3, while that of relative intensity decreases by a factor 
of about 4 from Fe X to Fe XV. These oscillations may be caused by the leakage of the photospheric 
p-modes through the chromosphere and transition region into the corona, which has been suggested 
as the source for intensity oscillations previously observed by TRACE. The temperature dependence 
of the oscillation amplitudes can be explained by damping of the waves traveling along the loop with 
multithread structure near the footpoint. Thus, this property may have potential value for coronal 
seismology in diagnostic of temperature structure in a coronal loop. 

Subject headings: Sun: atmosphere — Sun: corona — Sun: oscillations — waves — Sun: UV radiation 



1. INTRODUCTION 

A variety of propagating and standing MHD waves 
(e.g., slow mode, Alfven, and fast mode) have been ob- 
served in the Sun's outer atmosphere. They are mainly in 
coronal loops, but also in other structures such as coronal 
pluni es and proni inence s (see reviews by AschwandenJ 
2QQ3: IWajigl^QOl [2QQ5I: iNakariakov fc Verwichte h2005[ 
De Moortel I l2005l : iBaneriee etldTI I2007D . The detec- 



Similar intensity disturbances propagating along 
active-region loops with speeds of 122±44 km s~^ and 
intensity variations of 4.1%ibl.5% were observed in Tran- 
sition Region and Cor onal Explorer (TRACE) Fe IX/X 
171 A data ([Nightingal e et al. 1ll999l : 
1 20001: iRobbrecht et a ll l200ll: 



■ iDe Moortel et al. . 

iMcEwan fc de Moortel I l2006l) and EIT/SOHO 



De Moortel et al. 



tion of MHD waves in the solar corona is crucial to de 
termine the presence and relevance of wave-based heat- 
ing mechanisms since they are an obvious candidate to 
transport energy from the solar surface into the solar 
atmosphere. Such observations may also be used to 
improve existing estimates of coronal properties, both 
from direct measurements and indirect methods such 
as coronal seismo l ogy (s e e, e.g., [ Robert s et al. I 
Nakariakov et al.~l Il999l : INakariak ov fc Ofman 



Roberts fc Nakariakov 
20051 )! 



1984 : 



2001 



2003r INakariakov fc Verwichtel 



The first detection of propagating slow magnetoacous- 
tic waves was made in Ultraviolet Coronagraph Spec- 
trometer (UVCS) on board the Solar and Heliospheric 
Observatory (SOH O) observations of coronal h oles, high 
above the limb bv lOfman et al. I (|1997l . l2000ah . Similar 
compressive disturbances, with amplitudes of the order 
of 10—20%, and per iods of 10—15 min, were de tected 
in polar plumes by iDeForest fc Gurman I (|l998l ) with 
Extreme-ultraviolet Ima^in gs; Telescope (EIT)/SOHO. 
iQfman et alTI (|l999l . l2000b ^ identified the observed com- 
pressive disturbances as propagating, slow magnetoa- 
coustic waves, damped by compressive viscosity. 

^ Visiting Associate Professor, Tel Aviv University, Israel 



Fe XII 195 A data (|Berghmans fc Clette lfT999l ). These 
disturbances were also interpreted as slow magnetoacous- 
tic waves (Nakariakov et al. 2000). De Moorte fet al. I 
(|2002af ) found a distinct difference in dominant periods 
of waves in loops situated above sunspots (172±32 
s) and those above plage regions (321 ±74 s). This 
difference suggests that these waves probably origi- 
nate from the underlying oscillations, i.e., the 3 min 
chromospheric/transition-region oscillations in sunspots 
and the 5 min solar global oscillations (p-modes). 

Many authors have found evidence for the existence 
of 3 min sunspot oscillations propagating through the 
chrom osphere and tra nsition region into the lower corona 
(e.g. [B rvnildsen et al. Ill99 9a.b: Brvnildse et al. I l2002l : 
iFludra 2001: O ^Shea et al. .2002 : Rendtel et al. II2003D . 
The amplitude of the oscillations is found to reach a peak 
in the transition region lines then decrease with increas- 
ing tempe rature. These results are confirmed by a recent 
study of M arsh fc Walsh I ((2006') who clearly showed that 
the 3 min umbral transition region oscillations are di- 
rectly connected to the 3 min wave propagation along 
the TRACE loops. 

The p-mode oscillations are normally evanescent be- 
cause their periods are well above the cutoff period 
in the upper photosphere and low chromosphere as 
well as transition region in non-magnetic solar atmo- 
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sphere (e.g., lErdelvi et aDl2QQ7D. 

tic gravity waves, however. iBel & 



TABLE 1 
Emission lines observed with EIS. 



Fig. 1. — Left panel: The intensity map in the Fe XII 195.12 A hne from EIS. The raster observation was taken from 10:42 to 11:52 UT 
on 2007 Feburary 2. Right panel: The Doppler velocity measurements. The red color represents the red shift and the blue color the blue 
shift with a scale range from —20 km s"-*^ to +20 km s"-*^. The vertical line in each plot shows the position of the EIS 1 slit during the 
sit-and-stare observation, which was taken from 12:39 to 23:12 UT. The short horizontal lines on the slit mark the range of rows showing 
oscillations. The coordinate system for the figure corresponds to the end time of the raster observation. 

For m a ^netoa cous- 
Lerov I (|l977r ) pre- 
dicted that in regions of low-/3 plasma, the highly in- 
clined magnetic field can channel the low-frequency 
waves from the photosphere into the overlying coro- 
nal atmosphere due to the increase of the cutoff pe- 
riod. Numerical simulations based on this theory have 
demonstrated the appearance of 5 min waves in chro- 
mosph eric spicules and in coronal lo ops near active re- 
gions (|De Pontieu et al. I [2QQ3 . l2QQ5f ). The presence of 
this effect in the lower chromosphere was confirmed re- 
cently by ma ny authors with the TRACE 1700 and 
1600 A data (Mcintosh 



Jefferies "2006': " Jefferies e t al. I 
120061 : [Hoomfield et al. 2006; Vecchio et al. 2007 ). In 
particular, Fontenla et al. I ((1993a) and Jefferies e t al. I 
|200£) pointed out the importance of the observed low- 
frequency (<5 mHz) waves which may provide a signif- 
icant source of energy for heating solar chromosphere. 
Moreover, the detection of these waves propagating from 
the low atmosphere into the corona is also needed for our 
understanding of the energy balance in the outer solar 
atmospher e. Unfortunately, suc h observations are very 
few. Only 'Marsh et al. ' (12003) report the observation 
of 5 min propagating oscillations simultaneously at chro- 
mospheric, transition region and coronal temperatures. 

The EUV Imaging Spectrometer (EIS) onboard Hin- 
ode can simultaneously capture many emission lines 
from transition region to coronal temperatures, provid- 
ing us a good opportunity to study the temperature- 
dependent behavior of the oscillations and the propaga- 
tion of waves in the solar atmosphere (Mariska et al. I 
I2008D and is hig:hly valuable f o r coronal seismol- 
ogy ([yan Doorsselaere et al. I l2008l : lErdelyi fc Taroyan I 
l2008l ). In this study we report the first Hinode/EIS ob- 
servation of the 5 min upwardly propagating slow magne- 
toacoustic waves simultaneously at transition region and 
coronal temperatures in a plage region. The oscillations 
show up in both Doppler shift and line intensity with 



Ion 


Wavelength (A) L 


Og Tmax(K) 


He II 


256.32 


4.90 


Fe X 


184.54 


5.98 


Si X 


258.37/261.04 


6.13 


Fe XII 


195.12 


6.13 


Fe XIII 


202.04 


6.20 


Fe XIV 


264.78/274.20 


6.28 


Fe XV 


284.16 


6.32 


Fe XVI 


262.98 


6.43 



an amplitude dependent on the temperature. Section 2 
describes the observations. Two examples of oscillation 
packets are analyzed in Sect. 3 and the properties of 
their temperature dependence are presented in Sect. 4. 
Interpretation and discussion are given in Sect. 5, and 
finally our conclusions in Sect. 6. 



2. 



OBSERVATIONS 
2.1. Data 

EIS 



An overall description of EIS is available in 
ICulhane et al. I (l2007f). and th e Hinode mission is de- 
scribed bv lKosugT et al. I (|2007h . EIS has both imagmg 
(40'' and 266'' slots) and spectroscopic (1" and 2" slits) 
capabilities, in the wavelength range of 170-210 A and 
250-290 A with high spectral (0.0223A) and spatial res- 
olution (1" pixel"^). Its spectroscopic mode can operate 
in a rastering mode (repeated exposures while scanning 
over the observation target) or a sit-and-stare mode (re- 
peated exposures at the same spatial location). 

The observations analyzed in this study cover the cen- 
tral portion of NOAA active region 10940 and were ob- 
tained on 2007 February 2, when it was located close 
to the disk center. An EIS spectroheliogram was taken 
from 10:42 UT to 11:5 2 UT with the 1" slit and covering 
a 256" X 256" region (|Warren et al. II2008D . The expo- 
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sure times were 15 s. This data set contained 20 spectral 
windows. 

A sit-and-stare observation within the region began at 
12:39 UT and consisted of 1200 exposures with the 
slit, each with an exposure time of 30 s. Each exposure 
covered 20 data windows on the EIS detectors. This pa- 
per presents results for the lines in ten of those windows. 
Each window was 24 spectral pixels wide and covered 
a height of 400^' in the north-south direction. Table [1] 
lists the emission lines included in this study and their 
temperatures of formation. 

The raw data were processed by the standard Solar- 
Soft routine eis-prep to remove detector bias and dark 
current, hot pixels, and cosmic rays, and calibrated using 
the prelaunch absolute calibration. The EIS slit tilt and 
orbital variation in the line centroids were also removed 
from the data. The emission lines in each spectral win- 
dow were then fitted with Gaussian profiles, providing 
the total intensity, the Doppler shift, and the line width. 
The data in the short-wavelength detector were shifted 
downward by 18.5 pixels in the y-direction to correct for 
the offset between the two detectors. 

Figure [1] shows the intensity and Doppler shift maps 
of the active region covered by the EIS spectroheliogram 
in the Fe XII line. A string of brightenings in a circular 
shape are dominated by the redshifted flow. The location 
of the EIS slit for the sit-and-stare observation is across 
two loop systems. The large one at about y=0'' — 100'^ 
is most clearly seen in the Fe XV and Fe XVI lines (see 
IWarren e t al. 11200 8'). indicating that it consists of hot 
loops with the temperature of 2— 3 MK. The patches of 
brightening seen in Fe XII appear to locate at the foot- 
point region of these hot loops. 

The small one at the region from y=—60" to y=—50" is 
most clearly seen in the Fe XI— Fe XIV lines, indicating 
a temperature of 1—2 MK. The sit-and-stare observa- 
tion, made at the footpoint of this small loop, reveals 
oscillations in both Doppler shift and intensity in many 
emission lines. We will present the detailed analysis of 
the oscillation in Sects. [3] and [H Data from the GOES 
X-ray monitors show no flare-like fluctuations above the 
B level during the EIS observation. 

2.2. Effect of instrumental jitter 

Hinode is known to have instrumental jitter in both x 
and y directions that can have a range of up to 4'' during 
the observation. We first examine the jitter effect during 
the sit-and-stare observation. The left panels of Figure [2] 
show time series of the intensity and measured Doppler 
shift for the Fe XII line along the EIS slit. We found that 
three evident "dip" -like features occurred at about 19:23, 
21:02, and 22:40 UT. Since such "dip" features are seen 
simultaneously at three bright structures at the positions 
of y ^ —60'^ 40'^ and 70'\ we can determine that they 
were caused by the jitter of several pixels in the y di- 
rection. No high frequency, large amplitude jitter in the 
^/-direction are found in other times. The quasi-periodic 
Doppler shift oscillation was found to be associated with 
the small loop at ~ —60''. The oscillating structure 
has a width along the slit over more than 10^' and the 
oscillation is visible over the entire observation. These 
facts indicate that the detected oscillation could not be 
caused by jitter in the y direction. 

However, if there is a high gradient in brightness and 



Doppler shift across the slit, jitter in the x direction 
may cause artificial oscillation. To examine if this is the 
case, we model the sit-and-stare observation based on the 
EIS spectroheliogram and the drift of the XRT pointing 
which is obtained using the SolarSoft routine xrt-jitter. 
The XRT pointing is a good proxy for EIS pointing al- 
though not perfect. Since no XRT data are available 
before about 18:04 UT, we choose to model and analyze 
the sit-and-stare data set taken after 18:04 UT. The top 
panel of Figure [3] shows the displacements of the XRT 
pointing determined from the housekeeping data. Both 
the X- and ^/-displacements are mainly orbitally varying, 
with amplitudes up to 2" -3". The sit-and-stare obser- 
vation is modelled by extracting the slit slices from the 
EIS spectroheliogram taken at 10:42 UT (see Fig.[Tj) with 
the position of the slit varying with time, whose drifts 
in X- and ^/-directions are taken as the displacements of 
the XRT pointing. Since the intensity and Doppler shift 
distributions are assumed not to change with time, the 
fluctuations shown in the constructed time series (right 
panels in Fig. [2|) are only caused by the instrumental 
jitter. The evident "dip" -like features at about 19:23, 
21:02, and 22:40 UT which were caused by large dis- 
placements in the ^/-direction are consistent well with the 
EIS observation, confirming that the XRT pointing is a 
good proxy for the EIS. The orbital variations in bright- 
ness at y=—20" were also well reproduced. This feature 
was caused by the slowly- varying, orbit-related displace- 
ments in the x-direction which led the EIS slit repeatedly 
approaching a small brightening in the west. 

The middle and bottom panels of Figure [3] show com- 
parisons between the observed and modelled time profiles 
of the intensity and Doppler shift averaged at the loop 
of interest over 11 pixels from y=—b7" to —A7" (see the 
marked positions in Fig.[2|). Obviously, the variations in 
Doppler shift caused by the jitter are too small to ac- 
count for the observed oscillation. Except at the times 
when the large ^/-displacements occurred, there is no cor- 
relation in intensity variations between the modelled and 
observed time profiles. Therefore, we can safely exclude 
the possibility that the detected oscillation was caused by 
jitter of the EIS pointing. In addition, for the EIS sit- 
and-stare data to be analyzed in the following sections, 
the pointing drifts in the y direction have been removed 
based on the ^-displacements observed in XRT. 

3. ANALYSIS OF OSCILLATION PROPERTIES 

The quasi-periodic oscillation in Doppler shift detected 
at the footpoint of the small loop is clearly seen to ex- 
ist during the whole observation (see middle panel of 
Fig. [3]). The average root-mean-square (rms) amplitude 
of the time series with the subtracted background trend 
is about 0.7 km s~^. It is clear that this oscillation 
was not related to any flares or impulsive energy release 
events since the GOES X-ray flux was below the B-class 
level during the observation. We notice that the quasi- 
periodic time series is characterized by a train of oscilla- 
tions with a nearly constant period on the order of 5 min 
and an amplitude on the order of 1 km s~^, and these os- 
cillations are associated with intensity fluctuations with 
the same period. For example, two typical oscillation 
periods are 18:10-19:15 UT and 21:54-22:47 UT. 

In the following we analyze these two oscillations with 
the wavelet method. The details of the procedure are 
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Fig. 2. — Top left: Time series of the EIS sit-and-stare intensity in the Fe XII 195.12 A hne. Bottom left: The corresponding Doppler 
shift measurements. Right panels: The constructed time series of the Fe XII intensity and Doppler shift, based on the raster image 
and spacecraft drifts determined from the housekeeping data for XRT. The horizontal lines in each plot mark the range of rows showing 
oscillations. The vertical lines mark two time ranges for oscillations analysed (i.e., 18:10—19:15 UT and 21:54—22:47 UT). The Doppler 
shifts are plotted in a scale range from —20 km s"-*^ to +20 km s"-*^. 

assuming the background spectrum as white noise (with 
a flat Fourier spectrum). 



3.1. Oscillation 1 (18:10-19:15 UT) 

The left panel (a) in Figure 2] shows the evolution of the 
Doppler shift in the Fe XII line averaged over 11 pixels 
along the slit between y=—h7" and —A7". In practice, we 
first subtract the slowly- varying background trend from 
the time series (left panel (b) in Fig. 2]). The trend is 
constructed by using the moving average method with 
a characteristic smoothing time of 10 min. Then the 
wavelet spectrum and the global wavelet spectrum are 
constructed for the detrended time series (left panels (c) 
and (d) in Fig. [3|). The global wavelet spectrum is the 
average of the wavelet power over time at each oscilla- 
tion period. Due to the limited temporal resolution only 
periods more than 1 min are considered. The wavelet 
spectrum shows strong power at the period in a range of 
4—6 min over a duration of 6 periods. We measure the 
oscillation period as the value where the global wavelet 
power is peaked and the uncertainty as the half FWHM. 
The obtained period, Pyi, is 5.2±0.9 min and the am- 
plitude (defined as the square root of the peak global 
wavelet power), Ayi^ is 1.9 km s~^. 

The right panel (a) of Figure [4] shows the evolution of 
the Fe XII intensity. The intensity time series has been 
normalized by the mean value. A similar wavelet anal- 
ysis is applied to the detrended time series (right panel 
(b) in Fig. [4]). The wavelet spectrum shows that most of 
the power is concentrated within two period bands rang- 
ing in 4—6 min and 10—15 min (right panels (c) and (d) 
in Fig. [4]). For the short-period band, we measured the 
oscillation period, P/i=5.2±1.0 min and the relative am- 
plitude. All ^5.0%. The oscillation period is the same 
as that for the Doppler shift. The strong power is also 
seen over a time in agreement with that for the Doppler 
shift oscillation. The long-period band has an oscilla- 
tion period of 12.4±2.2 min and a relative amplitude of 
8.7%. No significant power is found at this period in the 




20:00 21:00 22:00 

Start Time (02-Feb-07 18:10:02) 



Fig. 3. — Top panel: Drifts of the spacecraft pointing in x 
and y directions, obtained from the housekeeping data for XRT. 
Middle panel: Average time profiles of the observed and modelled 
Doppler shifts for a strip shown in Fig. [2] Here positive values for 
the Doppler shift represent the blueshifted emission. Bottom panel: 
Average time profiles of the observed and modelled intensities. In 
the middle and bottom panels the thick horizontal lines mark two 
periods, over which the oscillations are analyzed with the wavelet 
method shown in Figs. |4] and [5] 



given by iTorrence fc C ompo^ (1998*) . For the convolu- 
tion of the time series the Morlet wavelet is chosen, and 
to establish whether the oscillations are real, a random- 
ization method is implemented which estimates the con- 
fidence level of the peaks in the wavelet spectrum by 
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a) Doppler shift time series a) Fe XII intensity tinne series 




18:20 18:30 18:40 18:50 19:00 19:10 0.0 2.5 5.0 18:20 18:30 18:40 18:50 19:00 19:10 0.000 0.004 0.008 

Start Time f02-Feb-07 18:10:341 Power (km/sY Start Time f02-Feb-07 18:10:34') Power 



Fig. 4. — Wavelet analysis for averaged Doppler shift and intensity time series in the Fe XII 195 A line from 18:10 to 19:15 UT. Left 
Panels: (a) The Doppler shift data (solid line) and the background trend (dotted line), (b) The detrended Doppler shift data, (c) The 
wavelet power spectrum. The dark color represents high power and the dotted contour encloses regions of greater than 99% confidence 
for a white-noise process. Cross-hatched regions on either end indicate the "cone of influence", where edge effects become important, (d) 
The global wavelet spectrum (solid line) and its 99% confldence level (dotted line). Right Panels: Same as the left panels but for averaged 
intensity time series. 



a) Doppler shift time series 
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b) Doppler shift time seris with subtraoted background 
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Fig. 5. — Wavelet analysis for averaged Doppler shift (left panels) 
21:54 to 22:47 UT. The annotations are same as in Fig. |4] 

wavelet spectrum of the Doppler shift. 

3.2. Oscillation 2 (21:54-22:47 UT) 

Oscillation 2 also showed up clearly in both the 
Doppler shift and intensity in the Fe XII line (Fig. [5]). 
The strong power for the Doppler shift is located at a 
period range of 3—7 min, with a slight shift from shorter 
to longer period with time (left panel (c) in Fig. [5|). We 
measured the oscillation period, P\/2=4.4±1.0 min, and 



a) Fe XII intensity time series 




22:00 22:10 22:20 22:30 22:40 



b) The detrended Fe XII intensity time series 
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Time 

c) Wavelet Power Spectrum d) Global Spectrum 




and intensity (right panels) time series in the Fe XII 195A line from 

the amplitude, Ay 2 =2.0 km s~^, from the global wavelet 
spectrum using the same method as applied to oscillation 
1. The wavelet spectrum for the intensity also shows 
two strong power bands as oscillation 1, one covering a 
shorter period range of 4—6 min and the other a longer 
period range of 13—20 min (right panel (c) in Fig. [5]). 
In contrast to oscillation 1, most of the power for the 
long-period band of this oscillation is within the cone of 
the influence where edge effects become important due to 
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Fig. 6. — Phase relationship between Doppler shift and intensity 
oscillations in the Fe XII line for two wave packets from 18:10 to 
19:15 UT {upper panel) and from 21:54 to 22:48 UT {bottom panel). 
In the upper panel, the dotted curve represents the reconstructed 
intensity time series by subtracting the power of the long-period 
band (8—16 min) from the wavelet spectrum. Here positive values 
for the Doppler shift represent the blueshifted emission. 

the finite lengtii of time series. For the short-period os- 
cillation, the period, P/2, is 4.8±0.9 min, and the relative 
amplitude, A/2, is 5.4%. The wavelet analyses for the in- 
tensity and the Doppler shift show that the short-period 
oscillation occurs in the same frequency range and dur- 
ing the same period. For the long-period oscillation in 
intensity, the period is measured to be 16.1±3.2 min and 
the relative amplitude is 6.2%. The measured physical 
parameters for oscillations 1 and 2 are listed in Tabled 

3.3. Phase Relationship 

We examine the phase relationship between the 
Doppler shift and intensity oscillations. For oscillation 
1, we reconstruct the intensity time series by removing 
the long-period (12 min) power band from the wavelet 
spectrum. The top panel of Figure [6] shows that the re- 
constructed intensity time series (dotted curve) is nearly 
in phase with the original one, and has a phase earlier 
by about 1/4-period than the Doppler shift oscillation. 
The cross-correlation gives a phase shift of 53°. The bot- 
tom panel of Figure [6] shows a good in-phase relationship 
between the Doppler shift and intensity oscillations for 
oscillations 2. The phase analysis gives the phase shift 
of 18°. Notice that for both oscillations the phase of 
the intensity oscillation is earlier than the Doppler shift 
oscillation. Examination of other cases show that the 
approximate in-phase relation is more common. 

According to linear MHD wave theory, intensity and 
Doppler shift oscillations are usually associated with a 
slow magnetoacoustic longitudinal wave in coronal loops. 
From the EIS spectroheliogram the oscillation is detected 
near the footpoint of a small loop so that longitudinal 
motions should have a line-of-sight component resulting 
in the observed Doppler shift. The phase relation of os- 
cillation 1 may indicate the presence of a standing wave 
or two oppositely propagating waves at that time. While 
the approximate in-phase relation for oscillation 2 may 
indicate that the oscillations are dominated by an up- 



wardly propagating wave, but probably overlaid with a 
weak downwardly propagating wave, which caused the 
small phase shift observed between intensity and Doppler 
shift oscillations. 

The propagation direction of the wave is determined 
based on the following linear wave theory. Given the axis 
of z is orientated towards the observer (i.e., the Doppler 
blueshift [upward motion] takes positive values), the dis- 
turbed velocity of an upwardly propagating linear wave 
can be described in the form. 



v{z^t) = V'sm{kcst — kz)^ 



(1) 



while for a downwardly propagating wave the disturbed 
velocity is. 



v{z, t) = V'sm{kcst + kz)^ 



(2) 



where V is the amplitude, k the wavenumber and Cs the 
sound speed. In the above we assumed uj = kcs with 
positive cj, /c, and Cg. The linearized continuity equation 
is 

where po is the background density (a constant), and p' 
the density perturbation. From Eqs. ([1]), ([2]) and ([3]) it 
follows that for an upwardly propagating wave 



p'{z,t) 



V 



posm{kcst — kz), 



(4) 



and for a downwardly propagating wave 

p {z, t) = - (-^j Posin(/cCst + kz). (5) 

The above equations indicate that the velocity and den- 
sity perturbations are in phase for the upwardly propa- 
gating wave, while in opposite phase for the downwardly 
propagating wave. Clearly, the in-phase relation for os- 
cillation 2 is consistent with the upwardly propagating 
wave. 

In contrast, the standing acoustic wave shows a 
quarter-period phase relation, and such examples have 
been observed by th e SUMER spectrometer on SOHO 
(jWang et al. I l2QQ3al lbl). In addition, we can exclude 
the possibility that the observed oscillations are caused 
by a fast sausage- mode wave, because the fast sausage 
mode is characterized by short periods on the order 
of se veral— tens seconds in coronal loops (jAschwandenl 
[200I . which are not consistent with the dominant pe- 
riod (5 min) of the observed oscillations. 

The fast kink mode is nearly incompres sible, i.e., with 
negligible density perturbation, however, ICooper et al. I 
(pQQ3l) demonstrate that intensity variations can be pro- 
duced by the kink modes, depending on the viewing an- 
gle. Since the intensity of the emission is proportional to 
the column depth of the loop, the variation of the LOS 
column depth due to the effect of projection causes the 
variation of the intensity. They find that the observed 
amplitude increases with the decreasing wave length and 
the increasing displacement amplitude of kink perturba- 
tions, and also depends on the angle between the LOS 
and the axis of the structure. However, the condition 
in our case is not in favor of such an effect because of 
the very long wavelength (inferred from the period of 5 
min), the small displacement amplitude (inferred from 
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Fig. 7. — Time series of EIS sit-and-stare intensity data (left panels) and the Doppler shift measurements (right panels) in six coronal 
hnes, Fe X (184.54 A), Fe XII (195.12 A), Fe XIII (202.04 A), Fe XIV (264.78 A), Fe XV (284.16 A), and Fe XVI (262.98 A), and a 
transition region hne, He II (256.32 A). Here the red color represents the redshift and the blue color the blueshift. The two horizontal lines 
in each panel mark the range of rows showing oscillations. The two vertical lines mark the time range for oscillation 2. The pointing drifts 
in y direction have been corrected in all time series. Note that the bad data points seen almost in all lines around 21:45 UT are caused by 
cosmic rays. 




^8:30 '8:40 '8:50 19:00 
Start Time (02-Feb-07 1 8 : 1 : 34) 



Start Time (02-Feb-07 21 :54: 10) 



Fig. 8. — Two selected time series of the averaged Doppler shifts (upper panels) and intensities (lower panels) in five coronal lines, Fe X 
(184.54 A), Fe XII (195.12 A), Fe XIII (202.04 A), Fe XIV (264.78 A), and Fe XV (284.16 A) for a strip marked in Fig. El The zero point 
of Doppler shifts in each emission line has been arbitrarily shifted in the y-direction in order to facilitate comparisons. The intensities 
in each emission line have been normalized to the mean averaged over the whole time series, and the light curves for the Fe X, Fe XII, 
Fe XIII, Fe XIV and Fe XV fines are shifted by -0.4, -0.2, 0, 0.2, and 0.4 in the y-direction, respectively. 



the Doppler shift amplitude) and the high inclination 
of the loop. Provided the LOS is co-planar with the 
plane of the kink oscillation and has an angle of 45° to 
the axis of the loop, we estimate that the intensity am- 
plitude produced by the kink oscillation should be less 
th an 1% with the measured parameters and the theory 
of ICooper et al. I (|2QQ3l ). which cannot account for the 
intensity amplitudes observed in our study. 

4. DEPENDENCE OF OSCILLATIONS ON THE 
TEMPERATURE 

4.1. Oscillations in Fe X—Fe XV 

We examine the temperature dependence of the oscil- 
lation in six coronal emission lines with formation tem- 
peratures ranging from 1.0 MK to 2.7 MK (see Table [T]). 



Figure [71 shows the evolution of intensity and Doppler 
shift for the oscillating loop in a part of time series. 
The background trend for Doppler shift time series has 
been subtracted at each position along the slit in order 
to show the oscillation more clearly. The loop is visible 
in all the six emission lines, but most clearly seen in the 
Fe XII— Fe XIV lines. The intensity oscillation is weak 
and not easily discerned, while the Doppler shift oscil- 
lation can be clearly seen in all lines except for Fe XVI. 
The quasi-periodic oscillation actually existed over the 
whole observation. At the time of oscillations 1 and 2, 
the oscillation appears to be more coherent along the slit 
(or across the loop) and more periodic compared to most 
of the other time. 
The top panels of Figure [8] show comparisons of the 
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22:10 22:20 22:30 22:40 22:50 23:00 
Start Time (02-Feb-07 22:04:43) 




22:10 22:20 22:30 22:40 22:50 23:00 
Start Time (02-Feb-07 22:04:43) 



Fig. 9. — Comparison between the Fe XII (195.12 A) and He II 
(256.32 A) time series in intensity (upper panel) and Doppler shift 
(bottom panel) from 22:05 to 23:08 UT. 

averaged Doppler shift profiles over 11 pixels along the 
slit between y=-h7" and -A7" in the Fe X-XV lines for 
oscillations 1 and 2. Strikingly, the plots clearly reveal 
a dependence of the oscillation amplitude on the plasma 
temperature, i.e., the amplitude decreases with increas- 
ing temperature. For oscillations 1 and 2, the maximum 
amplitudes in the Fe X line are about 2.5 and 3.1 km s~^, 
while those in the Fe XV line are 0.8 and 1.1 km s~^, re- 
spectively, indicating that the amplitude decreases by a 
factor of about 3. The oscillations seen at different tem- 
peratures are nearly in phase. Although the cross corre- 
lation shows that the phase of the oscillations in Fe X is 
slightly earlier than in Fe XV, the shift is measured to be 
within a half exposure time. The bottom panels of Fig- 
ure [8] show the evolution of the corresponding intensity 
profiles. For an easier comparison between the different 
emission lines, we have normalized the intensity for each 
line by the mean value of the time profile and shifted the 
light curve by a certain value along the y axis. Although 
the intensity oscillation is weak compared to the Doppler 
shift oscillation, the behavior of the temperature depen- 
dence of the amplitude is same as for the Doppler shift 
oscillation. For oscillation 1 three peaks can be seen in 
the Fe X and Fe XII lines, but are hardly discerned in the 
Fe XIV and Fe XV lines. For oscillation 2 the decreas- 
ing trend of the amplitude from Fe X to Fe XV is most 
evident for three peaks during 22:22-22:34 UT. 

The periods and the amplitudes for the Doppler shift 
and intensity oscillations in the five emission lines are 
quantitatively measured from the global wavelet spec- 
trum by applying the same method as used in Sect. 13. H 
and are listed in Table [2l The periods measured in the 
different lines are nearly same. For oscillation 1 the mean 
value of the periods for the Doppler shift in the five lines 
is 5.2±0.1 min, and for oscillation 2 the mean value of the 
periods is 4.6±0.2 min. The oscillation amplitude mea- 
sured for the Doppler shift from Fe X to Fe XV decreases 
by a factor of about 3. This result obtained with the 
wavelet method is consistent with that by directly mea- 
suring the maximum amplitudes mentioned above. For 



the intensity oscillations the measurement shows that the 
amplitude decreases by a factor of about 4. Note that 
the intensity oscillations in Fe XIII— Fe XV for oscillation 
1 and that in Fe XIII for oscillation 2 are so weak that the 
peak measured in the global wavelet spectrum is below 
the 99% confidence level. 

In addition, we can also examine the temperature de- 
pendence of the oscillation amplitude by measuring the 
average rms amplitudes. For the averaged time series 
of intensity and Doppler shift during 18:10— 23:10 UT. 
We first exclude bad data points at about 20:10 UT and 
21:45 UT which were caused by cosmic rays (see Fig. [7j), 
and then subtract the background trend which is taken 
as the 20-pixel smoothing average. Finally, the rms am- 
plitudes are calculated as the standard deviation for the 
mean of the detrended time series in five coronal lines 
and are listed in Table [2l We find that the oscillation 
amplitudes in Doppler shift and intensity decrease by a 
factor of 2.7 and 3.5, respectively, from Fe X to Fe XV. 
This result is in good agreement with that measured from 
oscillations 1 and 2 with the wavelet method. 

4.2. Oscillations in He II 

The 5 min quasi-periodic oscillations are also clearly 
detected in the transition-region line. He II 256.32 A, 
in both intensity and Doppler shift time series over the 
whole observation (see bottom panels in Fig.[7j). Interest- 
ingly, we did not find any correlation between the oscilla- 
tions observed in the Fe X— Fe XV lines (e.g., oscillations 
1 and 2) and those in the He II line, even considering the 
possible time delays. Figure [9] shows comparisons in the 
detrended relative intensity and Doppler shift between 
the Fe XII and He II lines. We have removed the con- 
tribution of a main blended line. Si X (256.37 A), from 
the He II line intensity by measuring the line intensity of 
Si X 261.04 A. Indeed, the blended emission from Si X 
can be safely ignored for the purpose of the oscillation 
study because it only contributes about 7% of the to- 
tal emission of He II in the data analyzed. Assuming 
that the Si X line has the same relative amplitude in in- 
tensity variation as the Fe XII line, we estimate that its 
contribution to the intensity variation in He II is below 
0.4%. The lack of correlation between the oscillations 
measured in Fe XII and He II can be clearly seen from 
Fig. [9] and the calculated correlation coefficients (see Ta- 
ble [3]) . The correlation coefficients between Fe XII and 
He II are about 10% in both intensity and Doppler shift, 
while those between Fe XII and Fe X are more than 50% 
in contrast. The similar result is obtained for a 5-hour 
period of data set. The above analysis indicates that the 
He II line detected the oscillation in transition region but 
not in the corona due to the blended emission. 

We analyze the oscillations seen in the He II line from 
about 19:30 to 20:30 UT with the wavelet analysis. Fig- 
ure [To] shows that two trains of oscillations can be iden- 
tified in the wavelet spectrum and their power distribu- 
tions in Doppler shift and intensity are similar (panels 
(c)). Each oscillation contains about three periods. We 
measure the oscillation period and amplitude from the 
global wavelet spectrum using the same method men- 
tioned above and list them in Table O We find that 
the period is in the range 4—6 min, and the Doppler 
shift and relative intensity amplitudes are comparable to 
those measured in the Fe X line. Figure [TT] shows that the 
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TABLE 2 

Measurements of the oscillations found in the EIS sit-and-stare data in six emission 

LINES^. 



Ion 


Pvi 




Pii 


An 


PV2 


Av2 


Pl2 


Ai2 


^(V) 






min 


km/s 


min 


% 


min 


km/s 


min 


% 


km/s 


% 


Fe X 


5.3±1.0 


3.1 


5.7±1.4 


6.8 


4.3±0.9 


3.2 


4.9±0.9 


6.5 


1.31 


3.8 


Fe XII 


5.2±0.9 


1.9 


5.2±1.0 


5.0 


4.4±1.0 


2.0 


4.8±0.9 


5.4 


0.72 


2.2 


Fe XIII 


5.1±1.1 


1.3 


(4.9±1.1) 


(2.2) 


4.9±1.2 


1.4 


(4.3±1.4) 


(2.6) 


0.69 


1.8 


Fe XIV 


5.4±1.7 


1.1 


(6.2±1.4) 


(3.1) 


4.6±1.1 


1.4 


4.4±1.0 


3.3 


0.63 


1.7 


Fe XV 


5.4±1.1 


1.2 


(4.5±0.9) 


(1.5) 


4.7±1.1 


0.9 


5.0±1.1 


2.0 


0.48 


1.1 


He 11^ 


5.7±1.0 


2.4 


6.1±1.1 


6.9 


4.4±1.0 


3.1 


4.8±0.9 


4.8 


1.07 


2.7 



^ Pvi and P/i are periods of the Doppler shift and intensity oscillations for oscillation 1 
(18:10-19:15 UT), respectively, while Pv2 and Pj2 for oscillation 2 (21:54-22:47 UT). Ayi and 
Aji are amplitudes of the Doppler shift and intensity oscillations for oscillation 1, respectively, while 
Av2 and A/2 for oscillation 2. The value in parentheses represents that the peak measured in the 
global wavelet spectrum is below the 99% confidence level. (j(V) and a{I) are the average rms 
amplitudes for the detrended time series in Doppler shift and relative intensity during 18:10— 23:10 
UT. 

^ Note that the two oscillations measured in the He II line are different from those measured in 
Fe X-Fe XV, which correspond to the time series from 19:30 to 19:52 UT and from 20:02 to 20:20 
UT. 




19:30 19:40 19:50 20:00 20:10 20:20 20:30 

Start Time (02-Feb-07 19:29:18) 

Fig. 11. — Phase relationship between Doppler shift and inten- 
sity oscillations in the He II line for time series from 19:29 to 20:32 
UT (upper panel). Positive values for the Doppler shift represent 
the blueshifted emission. 

Doppler shift and intensity oscillations are approximately 
in phase. Their phase shift measured with the cross cor- 
relation is 28°. Interestingly, the phase of the intensity 



TABLE 3 

Linear Pearson correlation coefficients between oscillations 
IN Fe XII AND Fe X lines and those in Fe XII and He II lines 

FOR INTENSITY (p/) AND DOPPLER SHIFT (pv) OVER A 1-HOUR AND A 
5-HOUR TIME SERIES. 



22:05 - 23.08 18:10 - 23:13 



Lines Fe XII/Fe X Fe XII/He II Fe XII/Fe X Fe XII/He II 

PI 0.52 0.09 0.57 0.04 

PV 0.61 0.12 0.53 0.16 



oscillation is also slightly earlier than that of the Doppler 
shift oscillation as found in Fe XII. The approximate in- 
phase relation between intensity and Doppler shift oscil- 
lations indicates the presence of an upward-propagating 
slow magnetoacoustic waves in transition region. The 
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TABLE 4 

Estimates of the wave energy flux for two oscillations 
measured in six emission lines^. 



Ion 






Fi 




F2 




km/s 


km/s 


10^ erg/s/cm^ 


km/s 


10^ erg/s/cm^ 


Fe X 


149 


5.1 


7.6±1.5 


4.8 


6.9±1.4 


Fe XII 


173 


4.3 


6.4±1.3 


4.7 


7.5ibl.5 


Fe XIII 


191 


2.1 


1.7±0.3 


2.5 


2.4ib0.5 


Fe XIV 


210 


3.3 


4.4±0.9 


3.5 


5.0±1.0 


Fe XV 


220 


1.6 


1.2±0.2 


2.2 


2.1±0.4 


He IP 


43 


^2.4 


^6.2 


^3.1 


^10. 



^ Cs is the sound speed corresponding to the hne formation tem- 
perature. V-^ and V2 cire the perturbed velocity amplitudes derived 
from the measured relative intensity amplitudes for oscillations 1 
and 2 (see Table 2). Fi and F2 are the wave energy flux calculated 
with Eq. 

^ Note that for the He II line the low limit of the perturbed ve- 
locity amplitudes is directly taken from the measured Doppler shift 
amplitude. 

small phase shift may imply that it is overlaid with a 
weak downwardly propagating wave as suggested for os- 
cillations seen in the coronal lines in Sect. 13.31 

5. DISCUSSIONS 

In this paper we have reported the simultaneous 
detections of 5 min quasi-periodic oscillations in the 
transition-region and coronal lines by Hinode/EIS in a 
corona loop rooted at plage. The oscillations show up in 
both Doppler shift and intensity throughout the whole 
observation, characterized by a series of trains with a 
more or less constant period lasting for 3—6 periods. 
These oscillation trains are more evident in Doppler shift 
than in intensity. The oscillations detected in five coronal 
lines in the range Fe X— Fe XV show a high correlation 
and exhibit a temperature dependence of the amplitude 
in both Doppler shift and intensity. The wavelet anal- 
yses show that the oscillation power is concentrated at 
the period in a range of 4—6 min. Both measurements 
with the wavelet analysis and the average rms method 
show that the oscillation amplitude decreases by a factor 
of about 3 in Doppler shift while decreasing by a factor 
of about 4 in intensity from Fe X to Fe XV. The oscilla- 
tions measured in the transition region line, He II, also 
have the dominant power at the period band of 4—6 min, 
and the amplitudes are comparable to those measured in 
the Fe X line. No correlation between the oscillations in 
Fe X— Fe XV and He II is found. The cross correlation 
shows that the phase of intensity oscillation is slightly 
earlier (by about 20° — 30°) than the Doppler shift oscilla- 
tion in Fe X— Fe XV and He II. The approximate in-phase 
relation indicates the presence of upwardly propagating 
slow magnetoacoustic waves in both the transition region 
and the corona near the footpoint of a loop. 

5.1. The source of waves 

From Eq. (|4|) we obtain p' / pQ = V /cg. We exam- 
ine this relation for oscillations 1 and 2 in Fe XII with 
the measurement obtained by the wavelet analysis, which 
gives the Doppler shift amplitude, V|| ^2.{) km s~^ and 

the relative intensity variation, /y/o ^5.2% (see Ta- 
ble [2]), where /q is the background intensity of the loop. 
Considering I' /Iq ~ 2p' / p^ and taking Cg = 170 km s~^ 
for the formation temperature of Fe XII, we deduce 



the expected value of the velocity amplitude, V ~4.5 
km s~^. The result of V|| = 0.43^ is consistent with the 

fact that the coronal loop, along which the slow- mode 
waves propagate, is highly inclined as seen in Fe XII (see 
Fig. [TJl). Assuming that is the LOS component of V\ we 
estimate that the inclination angle of the magnetic field to the 
vertical is about 65°. By applying the measured average rms 
amplitudes, almost the same value for the inclination angle 
is obtained. The high inclination of coronal fields in the loop 
may explain the absence of correlation between the oscilla- 
tions detected in Fe X— Fe XV and He II lines supposing the 
waves travel from the transition region to the corona along the 
loop, because the waves detected in the coronal lines and the 
transition region line should come from the different sources 
due to the inclination of the coronal fields. 

The oscillations we present here show many properties 
which are very similar to those of upward-propagating waves 
observed in coronal loops associated with plages u sing the 
TRACE 171 and 195 A bandpasses (e.g .Be Moortel et al. I 
2002a, tj; iKing et al.1 120031 : iMarsh et al. 1120031 1. First, these 
oscillations both have the dominant period of about 5 min 
and the amplitude of intensity variations of about 4%— 5% 
(in the Fe X and Fe XII lines for EIS). Second, they are both 
long-term existing and clearly nonflare-excited. Third, they 
both appear in the footpoints of the highly inclined coronal 
structures and are detected simultaneously at transition re- 
gion and coronal temperatures. These similarities indicate 
that the source of these waves is same, i.e., the leakage of 
the p-modes through the chromosphere and transition region 
into the corona. In our case the oscillations typically con- 
taining wave trains of 4—6 cycles with no apparent temporal 
damping are also consist ent with th e property of the photo- 
spheric p- modes. Be Pon tieu et al. I (p004) have shown that 
the inclination of magnetic flux tubes can dramatically in- 
crease tunneling and may even lead to direct propagation of 
the p-modes along inclined field lines. 

5.2. Interpretation for the temperature dependence 

Using the EIS data, we have detected the upwardly prop- 
agating slow magnetoacoustic waves simultaneously in five 
coronal lines with the formation temperature in the range 
'-^l— 2 MK. The waves in different lines show highly correla- 
tion (almost in phase). Particularly, the temperature depen- 
dence of the oscillation amplitude is revealed for the first time. 
The feature that the amplitude decreases with the increasing 
temperature is in good agreement with the observations of 
3 min sunspot oscillations which show that the amplitude 
peaks in the transition-region lines then decre ases with in- 
creasing temp eratur e jF ludra 2001; Brvnildse et al. I l2002l : 
lO'Shea et aLll2002l : iMarsh k Walsh . .2006) . The local maxi- 
mum amplitude of waves at the transition re gion can be ex- 
plain ed by strong stratification driven growth ( Erdelvi et al. I 
I2OO7I ). while the decreasing of the amplitude at increasing 
temperatures may be relevant to dissipation of slow mag- 
netoacoustic waves in the corona. The propagating distur- 
bances observed in the TRACE EUV images were found 
to be damped very quickly and typic ally only detected 
in the first 3-23 Mm along the loop (|Be Moortel et al. I 
2002bll . Numerical simulation s of slow MHB waves by 
Be Moortel HoodI (|2003l . 120041 1 showed that a combination 
of thermal conduction and area divergence yielded detection 
lengths that are in good agreement with observed values. 
Klimchuk et al. (2004) further developed the model consid- 
ering a non- isothermal loop, and found again that thermal 
conduction plays an important role in quickly damping of 
the waves. Here three possibilities are proposed to account 
for the temperature dependence of the amplitude for the ob- 
served propagating waves. 

(1) Assuming that the slow mode waves propagate up- 
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wardly in a non- isothermal coronal loop, then the smaller 
amplitude detected at higher temperature lines is because 
the lower temperature line forms in a height lower than the 
higher temperature line and the waves undergo a quick damp- 
ing during the propagation from the low level to the high 
level. With this interpretation time delays are expected to 
exist between the oscillations detected in the lines of differ- 
ent temperatures. Cross-correlation analyses for oscillations 
1 and 2 show that the phase of Doppler shift time series for 
Fe XV is slightly later than that for Fe X, but the time de- 
lay is less than 1 exposu re (i.e., <3Q s). By using t he results 
of the sunspot model of lObridko Staude I (|1988l ). we esti- 
mate the height difference of ~6400 km in vertical direction 
for the Fe X and Fe XV lines. This value is on the same 
order as estimated from observed hei ghts of EUV lines and 
the coronal model for the active Sun ( ISimon Noves I [197^ : 
ISimon et al. I ll974). Taking the sound speed of 200 km s~^ 
and an inclination angle of ^-^60° of the magnetic field to the 
vertical, we estimate the expected time delay to be about 60 
s, which is inconsistent with the observed. On the other hand, 
the inclination of ~60° implies a projected distance of more 
than 15" on the disk for the part of the loop with a verti- 
cal height of 6400 km. Since the slit appears to sit across 
the loop (see Fig.lTJi), this means that the emissions detected 
with the l" slit in Fe X and Fe XV lines should not come 
from the same loop. Thus, this scenario should predict no 
correlation between oscillations detected in the lines of differ- 
ent temperatures (at different heights), which contradicts to 
the observation. In the following we propose the other two 
alternative scenarios based on the loop with multiple threads 
of different temperatures. 

(2) With simultaneous TRACE observations in 171 and 
195 A bandpasses, iKing et al. I ((2003) have revealed that the 
slow mode waves propagate outwards along coronal loops of 
multiple-temperature structure. They show that the corre- 
lation between time series of disturbances observed in the 
different bandpasses has a tendency to decrease with dis- 
tance along the structure. They suggest that the initially 
high correlation may infer the same mechanism for genera- 
tion of the disturbances observed at different temperatures, 
while the decreasing correlation along the loop may be due 
to phase mixing of the waves. We may explain the temper- 
ature dependence of the oscillation amplitude with a similar 
picture, since the footpoint of the loop analyzed in this study 
is clearly not isothermal as revealed by the EM loci curves 
for Fe X— Fe XVI (not shown). Provided the waves have the 
same amplitudes at the base of the corona and propagate 
along parallel threads of different temperatures, we expect 
the wave of smaller amplitudes detected in higher tempera- 
ture lines because the dissipation of slow MHD waves is higher 
at the hotter plasmas by thermal conduction and compres- 
sive viscosity, which have been interpreted as the dominate 
damping m echanism in coronal loops by many theoretical 
studies (e.g. [ Nakariako v et al. ' f2 000l: lOfman Wang |[200l : 
iDe Moortel L Hood 2004: Klimchuk et al. II2004D . The in- 
terpretation in this picture, however, still needs to explain 
the high correlation between the oscillations detected in the 
different coronal lines, because time delays are expected if 
the waves travel at different speeds in different threads of the 
loop having different temperatures. One explanation could 
be that the time delay is too short to resolve with the 30 
s cadence. For example, if assuming that the waves travel 
over a distance of 7000 km (about half size of the footpoint 
brightening) in two threads with the temperature equal to 
the formation temperatures of Fe X and Fe XV, respectively, 
the expected time delay is only 15 s, less than the exposure 
time. 

(3) The third picture is proposed specifically to interpret 
the high correlation between the oscillations observed in coro- 



nal lines of different formation temperatures supposing that 
the absence of time delays is not due to the limited temporal 
resolution of the observation. We assume that the propa- 
gating waves are only present in an isothermal coronal loop 
with the plasma temperature of '-^l MK, which is overlaid by 
the hotter (>1 MK) loops with no waves propagating inside. 
The overlying relatively hotter loops are shown in the spec- 
troheliogram in the Fe XV line (not shown). Then the high 
correlation may be explained by the effect of the line response 
function (emissivity as a function of the temperature) since 
the modulated emissions of the different lines come from the 
same plasma disturbed by the waves. For the lines with higher 
formation temperatures, the emission contributed from the 1 
MK oscillatory plasma becomes less while those from the rela- 
tively hotter non-oscillatory plasmas become more dominant, 
therefore, the amplitude of the detected oscillations in both 
intensity and Doppler shift tend to decrease with increasing 
temperature of the emission lines. However, since both the 
Fe X and Fe XV lines form in a narrow temperature range, 
more exactly, the contribution of the emission to the Fe XV 
line from a plasma at 1 MK is more than 3 orders of magni- 
tude lower than its peak emission, the presence of waves only 
in the 1 MK plasma hardly explains the oscillation amplitude 
of 1%— 2% measured for Fe XV. In addition, this picture is 
lack of the theoretical basis for the assumption that the slow 
mode waves are only allowed to propagate in the 1 MK cool 
loops. 

Therefore, the second scenario provides the best interpre- 
tation of the temperature dependence of the oscillation am- 
plitude based on the present observations. Higher cadence (< 
15 s) observations are suggested in the future to check the re- 
Uability and accuracy of time delays between the oscillations 
detected in coronal lines of different temperatures, which will 
help us confirm this scenario. 

Based on linear wave theory we can explain why the ampli- 
tude of relative intensity decreases faster than the amplitude 
of Doppler shift with increasing temperature for the oscilla- 
tions we observed. From the relation of I' /2Iq ^ V /cs it 
follows that 



/y/o 

v 



(xT-'/^ (6) 

where T is the plasma temperature of the loop. If assuming 
Ti = l MK and T2=2 MK, respectively, as the temperature of 
the Fe X and Fe XV lines, we estimate the following ratio 
theoretically 



Rthe 



{I'/I,)tJ{I'/Io: 



T2 



T-2 



-1/2 



1.4. 



(7) 



This estimate is independent on scenarios proposed for ex- 
plaining the temperature dependence of the oscillation since 
Eq. (|7| is derived from the continuity equation. If assum- 
ing that the inclination angle of the field line for oscillations 
detected at T\ and at T2 is same, we estimate this ratio 



{I'/h)Tj{I'lh)T, {I'/Io)tJ{I'/Io)t, 



1.3, 



(8) 

with the observational measurements of 
{I'/Io)tJ{I'/Io)t, « 4 and {V{^)tJ {V(^)t, « 3. A 
good agreement between Robs and Rtheo provides fur- 
ther support to our interpretation of the observed 5 min 
oscillations in terms of a slow magnetoacoustic wave. 

For the oscillation in He II, we find that the amplitudes 
of relative intensity and Doppler shift measured with both 
the wavelet method and the average rms method are incon- 
sistent with the relation of I' /2Iq ^ P / Po — V /cg. For 
example, taking Cs=43 km s~^ we derive the perturbed veloc- 
ities are 1.5 and 1.0 km s~^ for oscillations 1 and 2, respec- 
tively, from the measured relative intensities (see Table [2]) , 
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which are evidently smaher than the measured hne-of-sight 
velocity amplitudes (2.4 and 3.1 km s~^). This inconsistency 
may imply that the measurements of relative intensity am- 
plitudes are underestimated, which could be relevant to the 
complexity of the He II line formation (e.g. iFontenla et al. I 
[l993b). Instead, we may estimate the low limit of the true 
amplitude for relative intensity from the measurements of 
Doppler shift amplitudes. From I' / Iq ^ 2yj|/cs, we ob- 
tain I' / Iq ^11% and 14% for oscillations 1 and 2, respec- 
tively. The estimated amplitudes are in g ood agreement 
with that measured by iMarsh et al. I (|2003l ) for the 5 min 
oscillation at the footpoint of fan- like TRACE loops. They 
find amplitudes of 12.4±2.1% for th e transition reg i on line , 
O V, observed with SOHO/CDS. Brvnil dsen eraLl ^mk) 
present SOHO/SUMER observations above a sunspot region 
and find intensity amplitudes of 11% and Doppler velocity 
amplitudes of 2.7 km s~^ for O V. Therefore, our measure- 
ments of Doppler shift amplitudes and the derived intensity 
amplitudes for He II are consistent with the amplitudes for 
O V found by other studies. Thus, our observations show that 
the relative intensity amplitude of 5 min oscillations decreases 
from the transition re^on to the c orona. This agre e s wit h 
the observations of l Fludral (l2001f): [ B rvnilds e et al. I (|2002l ): 
lO'Shea et afl (|2002r ) : IMarsh Wals h (2006) who find that 
the oscillation amplitude above the umbra reaches a maxi- 
mum for emission lines formed close to (1— 2)xl0^ K, and 
decreases for higher temperatures. 

5.3. Estimates of wave energy flux 

We can estimate the energy flux for the propagating waves 
measured in Fe X— Fe XV and He II lines in the coronal loop 
by, 

F = ^poiV'fcs, (9) 

where the sound speed is taken as Cs — 152T^^^ km s~^ 
with the plasma temperature, T, in units of MK, and V' — 
{l/2)cs{I' /Iq). For the oscillations observed in Fe X-Fe XV, 
we estimate V' from the relative intensity amplitudes mea- 
sured with the wavelet method, while for the oscillations in 
He II, we directly take the measured Doppler shift amplitudes 
as the low limit of V' for the reason discussed above. We use 
Si X A258.37/A261.04 and Fe XIV A264.79/A274.20 density 
sensitive ratios to diagnose the loop density, where the line 
ratio data are calculated with SSW/CHIANTI v5.2.1. For 
the data observed from 18:10 to 23:10 UT, we obtain the 
mean electron density as (2.0=b0.4) x 10^ cm~^ with the Si X 
line ratio, and (1.7±0.2) x 10^ cm"^ with the Fe XIV line ra- 
tio. We find that the values estimated from Si X and Fe XIV 
agree well. Note that in measurements of the Fe XIV A274.20 
intensity, the emission of a blended line Si VII A274.18 has 
been removed by considering the density-insensitive line ratio, 
A274.18/A275.35. We find that the blended emission from the 
Si VII A274.18 indeed can be ignored in this study because 
it only contributes about 3% to the Fe XIV A274.20. Using 
y9o=1.2mpiVe = (4±0.8)xlO~^^ g cm~^ (where iVe is taken as 
the electron number density measured with Si X line ratio, rrip 
the mass of proton, and a constant of 1.2 due to consideration 
of coronal He abundance) , we estimate the energy flux for os- 
cillations 1 and 2 in five coronal lines and list the values in Ta- 
ble |4l In estimates of the wave energy flux for He II, a typical 
value of mass density for the transition region, po=5xl0~^^ 
g cm~^, is taken. We find that the energy flux of waves for 
He II and Fe X is on the order of 10^ erg s~^cm~^, which 
decreases to the order of 10^ erg s~^cm~^ for Fe XV. Since 
the coronal radiative energy losses are ty pically on an o rder of 
10^-10'^ erg s ^cm ^ for active regions (|Aschwanden | [2004). 
the energy carried by the observed waves is too small to heat 
coronal loops. 



6. CONCLUSION 

In conclusion, the upwardly propagating slow magnetoa- 
coustic waves with periods of about 5 min have been detected 
in the transition region and coronal emission lines by Hin- 
ode/EIS at the footpoint of a coronal loop rooted at plage. 
The amplitude of the oscillations decreases with increasing 
temperatures. The temperature dependence of the amplitude 
observed in coronal lines can be explained by the waves trav- 
eling along a loop with multi-thermal temperature structure 
near the footpoint, and thus this feature may be valuable for 
coronal seismology to diagnose the property of multithreads 
in the loop. Many similarities between the waves observed by 
EIS in this study and the waves observed by TRACE in large 
fan-Uke loops suggest that the source of the waves is the same, 
i.e., a leakage of p- modes through the temperature-minimum 
region into the chromosphere and transition region reaching 
the corona. Although the energy carried by these waves is 
not enough to heat the corona, they may be important for 
the heating of the chromosphere, which have been found to 
be the sour ce for generation of th e periodic spicules in active 
regions (fDe Pontieu et al. ll2QQ? ). 
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